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Abstract 

Bose-Einstein correlations of charged kaons were measured near mid-rapidity in 
central Pb+Pb collisions at 158 A-GeV by the NA49 experiment at the CERN SPS. 
Source radii were extracted using the Yano-Koonin-Podgoretsky and Bertsch-Pratt 
parameterizations. The results are compared to published pion data. The measured 
m± dependence for kaons and pions is consistent with collective transverse expansion 
of the source and a freeze-out time of about 9.5 fm. 



1 Introduction 



Momentum correlations between pions have been widely used to study the 
space-time extent of the emitting source in nucleus-nucleus collisions. For iden- 
tical bosons the symmetry requirement of Bose-Einstein (BE) statistics results 
in an enhanced production of bosons with small momentum difference. This ef- 
fect was first observed by Goldhaber, Goldhaber, Lee and Pais (GGLP) [1] for 
pairs of like-charge pions at small opening angles. Kopylov and Podgoretsky 
[2] noticed the deep analogy with Hanbury-Brown and Twiss (HBT) space- 
time correlations [3] of the classical electromagnetic fields used in astronomy 
for interferometric measurement of star angular radii and developed the ba- 
sic methods of momentum interferometry in particle and nuclear collisions. 
This technique represents an important tool for investigating the underlying 
reaction mechanism as far as it influences the decoupling configuration. For 
example, an unusually long decoupling time was suggested as one of the sig- 
natures of quark-gluon plasma (QGP) formation [4,5]. 

In this paper, we present data from an experimental study of BE correlations 
between charged kaons emitted in central collisions of Pb nuclei at a beam 
energy of 158 A-GeV with a fixed Pb target. We compare the results with those 
from earlier analyses of pion data for the same reaction. It was argued [6] that 
correlations of kaons and pions might differ due to different contributions from 
resonance decays during their propagation toward freeze-out. A comparison 
could help to distinguish between a QGP deflagration and a conventional 
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Fig. 1. A schematic overview of the NA49 experimental apparatus showing the 
target T, the VTPCs within the superconducting magnets, the MTPCs and the 
TOF walls. Calorimeter VCAL provides the centrality trigger. 



production mechanism based on independent nucleon-nucleon interactions. 
However, rescattering between the final hadrons may smear out the differences 
and obscure what happened in the earlier stages of the collision. 

A probably more important aspect of kaons is associated with their higher 
mass. The particle source resulting from a collision between heavy nuclei is 
rapidly expanding in longitudinal and transverse directions. This manifests 
itself in a dependence of the extracted correlation parameters on the trans- 
verse mass m± of the particles, which form the pair [7]. Therefore, kaons and 
pions probe the situation in different kinematic regimes. Assuming that the 
decoupling surface coincides for both particle species, the kaons extend the 
m± range of the pion data and improve the sensitivity of comparisons with 
reaction models. 

Because of their much smaller abundance a study of kaon correlations requires 
higher event numbers and a significantly improved particle identification com- 
pared to the pion case. The only published results on kaon correlation in 
Pb+Pb collisions were obtained by experiment NA44 at the CERN SPS [8]. 
The NA49 spectrometer [9] includes a high-resolution time-of-flight (TOF) de- 
tector with acceptance at mid-rapidity which allows to select a relatively clean 
kaon sample. Data obtained from an analysis of approximately 1.7 million cen- 
tral as well as from 2 million semi-central events will be shown. The K + K + 
and K~K~ correlation functions are obtained in terms of the Bertsch-Pratt 
(BP) [4,5,10,11] and Yano-Koonin-Podgoretsky (YKP) [10,11,12] coordinate 
frames. They are corrected for Coulomb effects. The extracted parameters - 
correlation radii, freeze-out duration - will be compared with those obtained 
for pions [8,13,14,15,16] in the same central rapidity range. 
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Fig. 2. (a) Particle identification by combining TPC and TOF information. The 
circled region indicates the selected kaons. (b) TOF acceptance region in rapidity y 
and transverse momentum pt for single K + . Lines of constant momenta are shown 
by the dotted curves and the dashed line indicates the center-of-momentum rapidity 
ycM- (c) mj_-distribution for K + K + pairs. 



2 Experiment and Data 



The 158 A-GeV Pb ions from the CERN SPS were measured in position and 
direction by beam position detectors before hitting a Pb foil target of 224(336) 
mg/cm 2 thickness (1(2.5) % interaction probability) during the 1995/1996(2000) 
data taking periods. 

The trigger condition was derived from the signal detected in the veto calorime- 
ter (VCAL, see Fig. 1) which intercepts mainly projectile spectators 25 m 
downstream of the target. For the central data sample the 5% most central 
interactions were selected. In a geometrical picture this corresponds to impact 
parameters smaller than 3.5 fm. Data taken during different run periods con- 
tribute 380 000 (1995), 400 000 events (1996) and 911 000 events (2000). From 
the year 2000 run an additional sample corresponding to semi-central events 
(centrality 5% - 20%) was available and analysed. 

The NA49 spectrometer (described in detail in [9]) consists of two supercon- 
ducting magnets (total bending power of 7.5 Tm) with TPCs (time projection 
chambers) inside the magnetic field, two additional TPCs downstream of the 
magnets, and several TOF (time-of-flight) walls on both sides of the beam 
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(Fig. 1). The momentum vectors of charged particles are determined from 
their trajectories as measured in the TPCs and the position of the interaction 
point. The accuracy of the four-momentum difference^ q of a particle pair is 
of the order of 3.5 MeV in the momentum range of interest. 

For particle identification a combination of the specific energy loss dE/dx 
measured in the TPCs MTPC-R/L with the mass information obtained from 
the TOF walls TOF-R1/L1 was used. The method is illustrated by Fig. 2a for 
a typical momentum range. From Fig. 2a it is clear that no satisfactory kaon 
identification can be derived from either the specific ionization or the TOF 
signal alone. Therefore we restrict the phase space used in the present study to 
the acceptance of the TOF walls. Particles inside the region around the kaon 
peak (see Fig. 2a) are selected as kaons. The corresponding population in the 
rapidity (y) and transverse momentum (jpx) plane is displayed in Fig. 2b for 
the K + mesons; for K~ it is nearly identical. From fits of sums of Gaussian 
distributions for the different particle species the purity of the kaon sample 
is estimated to be about 93%, resulting in a purity of 86% for pairs. The 
distribution of the average transverse mass0 m± of accepted K + K + pairs is 
shown in Fig. 2c. Due to statistical limitations a subdivision into intervals 
was not possible. Because of slightly different detector efficiencies the weighted 
means of the correlation functions obtained from the 1995, the 1996 data and 
the 2000 data were used for the BE correlation analysis. Details of the analysis 
can be found in [17]. 



3 Bose-Einstein Correlation Analysis 



The correlation function C 2 (pi,p 2 ) is defined as the ratio of the probability of 
observing particles with four-momenta p± and p 2 simultaneously in one event 
divided by the probability for pairs of independent particles with the same 
single- particle phase space distribution. Experimentally, the raw correlation 
function is obtained as 



Cr{q)=N^ (1) 

where S(q) denotes the yield of particle pairs from the same event with the 
four-momentum difference q = p\ — p 2 . B(q) is the uncorrelated reference 
generated by the event-mixing method and the factor M is used to normalize 



1 Note that the units are chosen such that c = h = 1 

2 mj_ = \J m? K + Kj,, where Kt is the average transverse momentum of the kaons 
in the pair. 
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Fig. 3. 1-dimensional projections on of the 3-dimensional raw correlation functions 
Eq. 1 with qj < 36 MeV (j ^ i) in the BP parameterization. Top and bottom rows 
are for K + K + and for K~K~ , respectively. 

the correlation function to be unity at large relative momenta (200 - 2000 
MeV). 

Assuming an expanding Gaussian source in a representation given by Yano- 
Koonin-Podgoretsky [12] and Heinz [18] that is appropriate to explicitly take 
into account collective expansion, the correlation due to quantum statistics 
can be parameterized in terms of the components of q : 

C 2 (q±, q\\,qo) = 1 + \exp[-q\R\ 

-1yk{<1\\ - PYKqo?R\ ~ 1yk(Qo ~ PykQ^RI], (2) 

where q± = J (Ap x ) 2 + (Ap y ) 2 is the transverse momentum difference, q\\ = 
Ap z is the longitudinal momentum difference, qo = E\ — E 2 is the energy 
difference, A is the correlation intensity parameter and (3yk is the "Yano- 
Koonin velocity" which describes the source's longitudinal collective motion in 
the given rapidity range. The radii Ri, where i =_L, ||, denote the transverse, 
longitudinal and temporal extent of the coherence lengths, respectively. 

Although the YKP parameters have the most straightforward physical in- 
terpretation [19], this parameterization can be ill-defined in some kinematic 
regions [20]. Therefore, we use the Bertsch-Pratt parameterization in addition. 
For cylindrically symmetric systems and the midrapidity range it reads with 
vanishing cross terms [10] : 
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(3) 



where qi ong is identical with q\\ defined above, and (q ou t, Qside) are the com- 
ponents of the transverse 3-momentum difference in outward direction and 
perpendicular to it. As reference frame, a pairwise defined co- moving system 
would be appropriate. In view of the limited rapidity range considered here it 
is approximated by the so-called fixed local center-of-mass system, thus hav- 
ing the same frame for all particle pairs. Projections of the raw correlation 
function of K + K + and K~K~ pairs are shown in Fig. 3 in BP- variables. 

In order to extract the quantum-statistical effects, the counter-acting influence 
of the Coulomb repulsion as well as the dilution by misidentified pairs must be 
taken into account, while the effect of the strong interaction can be neglected 
[21]. This leads to the ansatz: 



where p = 0.86 denotes the purity of the pairs, A(q) represents the Coulomb 
repulsion factor and C 2 (g) is the correlation from quantum statistics param- 
eterized by Eqs. 2 or 3. Because of the considerable size of the source in the 
present reaction [22] the classical Gamov correction factor is not appropri- 
ate [23]. For pions the Coulomb correction can be approximately determined 
from the experimental correlation function for unlike-sign pairs [13], which 
is dominated by the Coulomb interaction and receives no contribution from 
BE correlation. In contrast, the strong interaction cannot be neglected [24] 
for unlike-sign kaons. Therefore, a parameterization of the Coulomb repulsion 
factor derived from a finite-size source Coulomb wave calculation [21] is used 
for the present analysis: 



where r* is the average distance between the emission points of the kaons in the 
pair center-of-mass system. Following [21] the value of r* can be calculated 
from the correlation radii Ri in Eqs. 2 and 3. As a consistency check the 
analogous procedure is applied below to the K + K~ correlations which are 
also measured in this experiment. 

The determination of the radius parameters Ri proceeds in an iterative way. 
Starting with estimated values for R iy the Coulomb repulsion factor A(q) 
(Eq. 5) is calculated and smeared with the measuring resolution of q. A set 
of Ri and A are then obtained from a 3-dimensional maximum likelihood fit 
of the r.h.s. of Eq. 4 to the raw correlation function Eq. 1. The procedure is 
iterated until the fit parameters reach stable values. 



(l-p)+p-C 2 (q)-A(q) 



(4) 



A(q)=A c (r*,q), 



(5) 
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Fig. 4. 1-dimensional projections on q^ of the Coulomb- and purity-corrected 3-di- 
mensional correlation functions Eq. 6 of like-sign kaons with qj < 36 MeV (j ^ i) in 
YKP (top) and BP (bottom) parameterizations. The solid lines in each histogram 
indicate the projections of the fitted 3-dimensional functions of Eqs. 2 and 3 with 
the parameters given in Table 1. 



Table 1 

Source parameters with statistical errors obtained from 3-dim fits of the central 
data ; the 1-dim projections are shown in Fig. 4. 



YKP 


A 


R\\[fm] 


R ± [fm] 


Rolfm] 


Pyk 


X 2 /NDF 


K+ 


0.74±0.06 


3.95±0.35 


4.34±0.30 


2.86±0.65 


0.01±0.02 


216/228 


K~ 


0.86±0.05 


4.52±0.45 


4.15±0.28 


3.00±0.52 


0.01±0.05 


270/328 


BP 


A 


Rlongifm) 


Rsideifm] 


Routifm] 




X 2 /NDF 


K+ 


0.81±0.06 


4.46±0.25 


3.58±0.40 


5.07±0.27 




439/446 


K~ 


0.89±0.07 


4.78±0.33 


4.55±0.31 


4.97±0.39 




499/537 



For the central collisions the final results for the YKP and BP parameter- 
izations are summarized in Table 1 with the fit quality given by \ 2 /NDF '. 
The quoted errors are only statistical. The systematic errors on the the fitted 
parameters due to uncertainty in Coulomb correction and strong interaction 
are estimated to be about 10% of each parameter. 

Using Eq. 4 and the fitted parameters one can calculate a Coulomb corrected 
correlation function C| orr from the directly measured raw correlation function 
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50 100 150 200 250 300 

QiJMeV] 

Fig. 5. K + K~ correlation function in central collision events plotted versus Qi nv . 
The curves show theoretical calculations (see text) with the effects of measuring 
resolution and misidentification applied. The dotted curve represents a pointlike 
source (Gamov factor). The finite source size Coulomb wave calculation is indicated 
by the dashed curve. The final result including also the effect of strong interactions 
is given by the full curve. 



C r 2 aw {q) -1+p 
p-A(q) 



(6) 



The result is shown in Fig. 4 as data points and compared to the fitted correla- 
tion functions which are depicted by the curves. One observes good agreement. 



In addition the kaons from about 2 million semi-central (5% - 20%) events 
were also analysed and the extracted parameters are listed in Table 2. Both 
the values of the radius parameters and of the correlation intensity are found 
to be smaller than for the central collision events, in qualitative agreement 
with the expected decrease of the radii with the decreasing centrality. 



Fig. 5 displays a plot of the K + K~ correlation function versus Q inv = J (Ap) 2 — 
for central collision events. The observed correlation is due to the Coulomb 
attraction as well as the strong interaction of the K + and K~~ [21]. The data 
points are compared to theoretical calculations for which the effects of momen- 
tum resolution and misidentification were taken into account. As demonstrated 
by the dotted curve, the Gamov factor (appropriate for a point like source) 
does not describe the data. The finite size Coulomb wave calculation using 
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Table 2 



Source parameters from the correlation functions of kaons from semi-central events. 



YKP 


A 


i?H [/m] 


R±ifm] 


Roifm] 


Pyk 


X 2 /NDF 




L).oo±U.Uo 


o o^in go 


z.9o±U.18 


9 A 1 _i_n TO 
0.4i±U. I Z 


n no_i_n no 
U.Uz±U.Uz 


o^^7Q /om n 
ZD M/zUiy 


K 


0.60±0.06 


4.07±0.51 


3.95±0.38 


2.03±0.82 


0.01±0.04 


1553/1574 


BP 


A 


Rlongifm] 


R s ide[fm] 


Routifm] 




X 2 /NDF 


K+ 


0.60±0.03 


3.35±0.35 


3.16±0.48 


3.40±0.26 




1529/2012 


K 


0.67±0.04 


3.55±0.50 


3.36±0.51 


3.95±0.37 




968/769 



the fitted radius parameters is shown by the dashed curve. It already pro- 
vides a much better representation, but still stays somewhat above the data. 
When also including the contribution from strong interactions (full curve) a 
good reproduction of the region Q inv < 200 MeV is obtained. This agreement 
can be taken as a confirmation of the Coulomb correction procedure used for 
like-charge kaon correlations. 



4 Discussion and Conclusion 

From Tables 1 and 2 one observes that the radius parameters for K + and K~ 
are similar with a trend to somewhat larger values for K - . This is consistent 
with the expected influence of Coulomb interaction with the positive charge 
of the reaction zone [22]. The correlation intensity parameter A stays below 
unity. Possible causes are e.g. a partially coherent source, a contribution to 
the observed kaons from decays of long-lived resonances, or a deviation from 
the Gaussian parameterization due to a mixture of sources with different radii 
[25] . The longitudinal source velocities (3yk are close to zero in agreement with 
the fact that our acceptance is restricted to midrapidity (see Fig. 2b). 

When comparing the radius parameters obtained from kaon correlations with 
those of pions one has to bear in mind that these parameters depend on the 
transverse momentum or transverse mass as a consequence of the collective ex- 
pansion of the source [7]. In Fig. 6 we show the earlier pion data [8,13,14,15,16] 
for the same reaction together with the kaon results from Table 1 as well as 
those from the CERN experiment NA44 [8]. The radii show mj a dependence 
with a ~ 0.5 for the longitudinal radii (as suggested by Makhlin and Sinyukov 
[26]) and a ~ 0.3 for the transverse radii. 

According to the relation : 
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Fig. 6. Comparison of radius parameters of kaons and pions as function of m±. Error 
bars show statistical and systematic errors added in quadrature. The dashed lines 
indicate the fit of the function Cm7° with different a, respectively. 




(7) 



based on a hydrodynamic expansion model [26], where 77 and T are the freeze- 
out time and temperature of the particles, respectively, the longitudinal radii 
can be plotted as a function of yT '/m± (Fig. 7a) using the known freeze-out 
temperature T ~ 120 MeV from [13]. From its gradient the freeze-out time 
is extracted as 9.5 ±1.5 fm, which is consistent with the previous estimate 
based on the pion correlation study in [13]. Furthermore, the relatively weak 
dependence of R± on m± can also be reproduced by hydrodynamic expan- 
sion models [10,13] with a transverse expansion velocity /3± around 0.55 (see 
Fig. 7b). 



In conclusion the kaon radius parameters are fully consistent with the pub- 
lished pion results and the hydrodynamic expansion model. The influence of 
different contributions of resonance formation and different rescattering cross 
sections seem to be either minor or similar in both cases. Pions and kaons thus 
seem to decouple simultaneously. 
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Fig. 7. (a) Rii as a function of y/T/m± according to Eq. 7. The solid line indicates 
the function with tj = 9.5 fm with ±1.5 fm indicated by the dashed lines, (b) R± 
as function of m±. The dashed line denotes a model calculation [10,13] with the 
transverse expansion velocity 0± = 0.55. 
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